A simple physical model of liquid-glass transition: intrinsic fluctuating interactions 
and random fields hidden in glass-forming liquids 
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We propose that glass-forming liquids are intrinsically under the influences of both fluctuating 
interactions and random fields well-known in the field of spin systems. This is due to the frustration 
between the isotropic and anisotropic parts of effective intermolecular interactions. Our model 
indicates the existence of two key temperatures relevant to glass transition, the density ordering 
point Tm and the Vogel-Fulcher temperature To. Between T^ and To, a system has features similar 
to the 'Griffiths phase', while below To it has those peculiar to the 'spin-glass phase'. This picture 
naturally and universally explains vitrification behavior from its strong to fragile limit. 

PACS numbers: 64.70.Pf, 61.43.-j, 64.60.Ht, 61.25.Em 



Although various features of the glass transition have 
been clarified recently there has so far been no 

simple physical description of glass transition covering 
its strong to its fragile limit. The theoretical approaches 
to this problem can be grouped into two types: (i) one fo- 
cuses on the slowing down of the dynamics on approach 
to a glass-transition temperature T g . The topological- 
confincmcnt effects on molecular motion have been ex- 
pressed in terms of a few different concepts, such as free 
volume Jj| , cooperative rearranging regions Q , and cage 
. This stream leads to the development of the mode- 
coupling theory (MCT) |Q. However, the physical mean- 
ing of assumptions hidden in mode-coupling approxima- 
tions is still not clear (ii) The other applies the 
knowledge of 'spin glass' P~|lO|, whose glassy behavior 
is much more deeply understood than that of structural 
glass, and 'frustrated systems' |pd| , p^| to the problem of 
glass transition. This type of approach is attractive in 
the sense that it has a potential to provide us with a 
universal physical picture of random frustrated systems. 
These two approaches are essentially different from each 
other in that the former presupposes disorder, while the 
latter puts more emphasis on the ordering phenomena. 
It should be noted that neither of the above approaches 
provides us with any clear answer about what physical 
parameter controls the fragility of liquids [jl]j^] . 

Before making any model of glass transition, we have 
to seriously consider a much more fundamental problem, 
namely, why some molecules crystallize without vitrifica- 
tion while others can easily form glasses without crystal- 
lization: Without considering 'crystallization', we cannot 
understand either a supercooled state, which is defined as 
a metastable state below a melting point, or vitrification 
phenomena, because the avoidance of crystallization is a 
prerequisite to them. Approach (i) regards glass tran- 
sition as a purely dynamic transition. Thus, it usually 
neglects even the fact that a glass-forming system has 
the ability of crystallization. Approach (ii), on the other 



hand, fails to give a clear molecular-level explanation on 
what is the physical origin of the similarity between spin 
glass and structural glass, although it recognizes the im- 
portance of topological frustration . Thus, the wide 
gap between them has not been filled yet. Hence, neither 
approaches can clearly explain why slow dynamics starts 
to appear when a liquid enters into a metastable super- 
cooled state upon cooling, namely, when a liquid is cooled 
below the melting point of its crystal. 

The main aim of this Letter is to clarify the physical 
origin of frustration in glass-forming liquids that appar- 
ently have no intrinsic quenched disorder and to estab- 
lish a direct connection between structural glass and spin 
glass on the basis of molecular-level consideration. In 
contrast to the common sense that density p is the only 
order parameter required for the physical description of 
liquids, we propose that the introduction is necessary of 
a new order parameter (bond parameter) representing a 
locally favored structure in liquids: This local bond for- 
mation causes 'fluctuations in the intermolecular inter- 
actions' and 'symmetry-breaking random fields' in much 
the same way as in spin systems ^,0- On the basis of 
the knowledge of spin systems [ flif , we propose a univer- 
sal picture of glass transition. To our knowledge, this 
is the first approach to the problem of the liquid-glass 
transition directly focusing on 'crystallization (or density 
ordering)'. 

First we focus on the effective attractive interaction po- 
tential between a molecule and its neighboring molecules. 
It is generally given by the form V(r, £1) = V(r) + 
AV(r, f2), where r is the distance from the center of mass 
of the molecule and 17 expresses the orientation. V rep- 
resents the isotropic part of the interaction and AV its 
anisotropic part. This anisotropy can lead to a locally 
favored structure made of a molecule and its neighbor- 
ing n molecules, which is not consistent with the crys- 
tallographic symmetry. Thus, there can exist competing 
interactions in any liquids. Even for spherical-particle 
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systems, it is known that an 'icosahedral structure' is 
locally favored, which is separated from alternative ar- 
rangements such as bcc and fee by high potential en- 
ergy barriers (~ nAV, n = 12). Its importance in the 
glass problem was widely recognized by many researchers 
[p|- |l2| , p"5| [l7| . This local bond ordering is incompatible 
with any long-range order, so that it plays the role of 
a random disordering field against crystallization: it fa- 
vors vitrification. This competition between density and 
bond ordering results from two conflicting requirements: 
(i) to minimize the distance between nearest-neighbor 
molecules and (ii) to maximize the number of surround- 
ing molecules. In strong liquids, on the other hand, the 
locally favored symmetry is mainly selected by specific 
anisotropic interactions between molecules 18 such as 
hydrogen and covalent bonding, which often lead to a 
'tetrahedral structure'. For a general physical descrip- 
tion of real liquids, thus, we need a new order parameter 
S to describe the presence of such locally preferred ar- 
rangement of liquid molecules: p and S are the minimal 
order parameters required for the physical description of 
the above complex features of many-body interactions. 
Although S should have a tensorial character that plays 
important roles in the selection of crystallographic sym- 
metry and its rotational dynamics, we can treat it as a 
scalar order parameter S when we consider the phase be- 
havior. The relevance of such an approximation is well 
established in the field of liquid crystals p9[ . Without 
loosing generality, therefore, the bond parameter S can 
be defined as the 'local number density of locally favored 
structures': S(r) = £j<5(r — fj), where f*j is the position 
vector of a locally favored structure (number i) and Ej is 
the sum about i over a unit volume. Note that the spatial 
distribution of locally favored structures (f\) is random 
due to the nature of bond formation [^o| . The spatially 
averaged value of S is given by S = So exp(ftnAV) , where 
(3 = 1/fceT (fee is the Boltzmann constant and T is the 
temperature), since a locally favored structure, which is 
stabilized by n bonds, is in a lower energy state than 
the other part of liquid by nAV. When the dependence 
of AV on is not consistent with the crystallographic 
symmetry, S decreases the local density, lowers the crys- 
tallization temperature and thus locally disturbs crystal- 
lization. Thus, each molecule itself internally has the 
cause of disorder and random fields against the density 
ordering. We believe that this frustration between p and 
S is the physical factor determining how easily molecules 
can form a glass without crystallization. Note that the 
strength of disorder effects increases with A V, and so the 
stronger glass suffers from the stronger spin-glass (SG) 
effects. 

The above consideration leads to the following two- 
order-parameter description of liquids: In contrast to 
the common sense that a liquid state can be described 
by a single order parameter, p, we need a new order pa- 
rameter S to express a tendency of local bond order- 



ing since in any liquids molecules locally favor a cer- 
tain symmetry, which is not necessarily consistent with 
a crystallographic symmetry favored by by density or- 
dering. This competition between the two order param- 
eters, p and S, causes frustration, which we believe is 
the major origin of vitrification. This physical picture 
can be naturally modeled by a Ginzburg-Landau-type 
model with the couplings between p and S, which rep- 
resent disorder effects of S on the ordering of p. Den- 
sity fluctuations Sp in the liquid phase indicating the 
instability toward the solid phase have a maximum at 
nonzero wavenumber qo, whose essential feature is well 
characterized by the following bare static structure fac- 
tor: S (q) = 2k B T X o(q) = k B T/[r + K(q - q ) 2 ]. The 
phenomenological free energy which predicts the above 
S (q) is given by |1 



iff, ■■ ■■ j drdr'Spirjxo 1 (r - r')5p(r') - | / ilrSp[ FY' 



+ | / drSp(f)\ 
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where r = 02 (T — T*) and T* is the temperature of 
the mean-field limit of stability of the liquid phase. This 
Hamiltonian implies that 'density order parameter favors 
spherical symmetry'. By including the coupling between 
Sp and S into the above standard theory of a liquid-solid 
transition, we obtain the following Hamiltonian that we 
believe is relevant to the glass transition: 

[3H GT = f3H P + J dr [- Cl Sp(r)S(f) - ^6p(r) 2 S(r)}. (2) 

In the above, the coupling between p and S is introduced 
through the coupling constants Cj. For the "negative" 
coupling < 0), the formation of active bonds (or S) 
leads to a decrease of local density and also to a decrease 
of the ordering temperature. For the "positive" coupling 
(cj > 0), on the other hand, the formation of active bonds 
(or S) leads to an increase in local density and to an 
increase of the ordering temperature, and so molecules 
should never form a glass and just crystallize. The type 
(sign) of coupling between p and S gives a simple crite- 
rion on whether molecules crystallize without vitrification 
or can easily form a glass. Hereafter we consider the case 
of Ci < 0, since vitrification can occur only for this case. 
The dynamics of Sp can be described by pel 



dSp{f,t) 
dt 



= r n v- 



5{PHgt) 
5(6p(r,t)) 



(3) 



where £ is the noise term and T is a bare kinetic co- 
efficient. Although S is not a quenched variable, the 
much slower dynamics of S than p guarantees the quasi- 
quenched nature of S: The lifetime of locally favored 
structures is longer than the characteristic time of den- 
sity fluctuations by at least a factor of R ~ exp((3nAV). 
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More importantly, recent theoretical studies indicate that 
even a frustrated random system without quenched dis- 
order has essentially the same features as that with 
quenched disorder ||. These facts allow us to regard 
S as a quenched variable as long as we consider disorder 
effects of S on density ordering. Thus, we treat S as a 
quenched variable hereafter. 

In relation to the above, we briefly consider the sta- 
bility of a metastable supercooled state. First we de- 
fine T^(0) = T* % and T^(S) respectively as the density- 
ordering temperature (a crystallization or melting tem- 
perature) without and with the effects of random dis- 
order (S is the disorder strength). The former is the 
melting temperature of a real defect-free crystal formed 
as a result of density ordering, while the latter is that 
of a hypothetical crystal with quenched disorder (or lo- 
cal bond ordering), which may not exist in reality. The 
melting temperature of this hypothetical crystal with de- 
fects, T* n (S), rapidly decreases with an increase in defect 
density, or S (see figure Q), which is known as 'dilution 
effects' in spin systems. For real crystallization, which 
requires the annealing of the random nature of S, on the 
other hand, a system has to overcome a large energy bar- 
rier corresponding to the break-up or deformation of all 
the active bonds in a critical nucleus of crystal, which 
costs the energy of the order of AE ~ nAVSv n , where 
v n is the volume of a critical nucleus. This naturally 
explains why a metastable supercooled state is so sta- 
ble in glass-forming liquids. Crystallization at can 
thus be kinetically avoided by this extra energetic bar- 
rier for nucleation, AE, for a sufficient cooling rate and 
thus vitrification can be induced, as described below. In 
the metastable branch of a supercooled liquid, therefore, 
the random distribution of locally favored structures is 
not altered [^0| by density ordering . Provided that crys- 
tallization is kinetically prohibited, we can regard a su- 
percooled state as the quasi-equilibrium thermodynamic 
state. Thus, we consider a problem of vitrification on the 
basis of this quasi-equilibrium assumption. 

Here we point out the similarity of the above Hamil- 
tonian Hqt [see equation (|J)] assuming a quenched na- 
ture of S* and that of a spin system under fluctuating 
interactions and random fields. The Ginzburg-Landau- 
type Hamiltonian under random transition temperatures 
<5r(r), Hqd, and that under random fields h(r), Hrf, 
are [f4| : 



Wqd = J df [^(V0) 2 + i(r - 5r(r))0 2 + ±g<j>% 

PH RF = J dr [i^(V0) 2 + ~r</> 2 + \g<f + h(r)4>], 

where 0(r) is the order parameter. Our Hamiltonian 
Hqt has both features of EIqd and Hrf at the same 
time: By setting (j){r) = Sp(f), 5r(f) — C2<5'(r), and 
h(r) = —ciS(r), Hqt can be directly correlated to Hqd 



and Hrf, after coarse graining the high q mode around 
go- Thus, the essential effects of random disorder on den- 
sity ordering in our system should be the same as those 
in spin systems The quasi- quenched nature of S 

allows us to consider the problem in terms of complex en- 
ergy landscape, without solving the dynamic equations. 

On the basis of the knowledge of spin glass IIH , we 
draw a simple physical picture of glass transition (see 
figure [l]). Above T* n , the system behaves just as an ordi- 
nary liquid and the structural relaxation time T a obeys 
a simple Arrhenius law, r a = exp(AG/fceT), where 
AG is the activation energy and an increasing function 
of AV, and the relaxation function $(i) is exponential: 
$(t) = $(0) exp(-t/r a ). Below T^, the saddle-point 
equation 5Hqt /8{8p) = 0, namely, —KV 2 Sp(r) + [r — 
C2S(r)]5p{f) + gSp(r) 3 — C2<S'(r) = 0, starts to have local 
minimum solutions around the high-density side of the 
liquid free-energy minimum pl[ . This means that there 
exist a macroscopic number of spatial 'islands' having a 
higher density than the liquid, below T* n . This situa- 
tion is similar to the spin glass where the phase space is 
factorized into a hierarchy of 'valleys', or pure states of 
local minima separated by macroscopic barriers. Thus, 
slow relaxations overcoming barriers separated by differ- 
ent valleys are expected. This phase existing between 
and T^(S), which is characterized by the existence 
of numerous metastable states separated by finite barri- 
ers, is similar to the "Griffiths phase" |22| , |23| ] known in 
spin systems. 

This 'Griffiths- like phase' can be characterized as fol- 
lows QIIQ. In the temperature interval T^(S) < 
T < T* l7 $(i) is described by the stretched exponential 
HQ: $(t) = $(0)exp[-(V'r Q ) /3ff ], instead of the usual 
exponential relaxation, as it should be in the ordinary 
liquid phase (the paramagnetic phase in spin systems) 
|P5| . In spin systems ||(| , it is claimed that the stretched 
exponential parameter (3k is the temperature-dependent 
exponent, as it is a finite value (< I) at T = T^S), and it 
increases monotonically up to (3k = I at T^. This is con- 
sistent with what has been observed in structural glasses 
HHH^] . Further, the strong coupling between 'islands' 
probably leads to the existence of the SG-like phase be- 
low Tsgi hi which the numerous disorder-dependent local 
minima are probably separated by macroscopic (or infi- 
nite) energy barriers [fblf . We assign this phase transition 
at Tsg from the Griffiths- like phase to the SG-like phase, 
to the Vogel-Fulcher temperature To, where the relax- 
ation time diverges due to the infinite barriers and the 
ergodic to non-ergodic transition takes place. However, 
this point is never reached in real experiments because 
the large barrier heights near To cause extremely slow 
relaxations even above Tq. Then, the glass-transition 
temperature T g can be defined as a temperature where 
the metastable 'islands' having sufficiently high energetic 
barriers do percolate. Thus, T* n > T g > T Q > T^(S). 

Figure [l] shows a phase diagram of glass transition on 
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the basis of the above physical picture. Here it should 
be noted that for the liquid-glass transition there can 
be two types of the origins of disorder whose tempera- 
ture dependencies are crucially different from each other: 
(i) anisotropic interactions that are not consistent with 
the crystallographic symmetry, in molecular glass form- 
ers ; and (ii) quenched disorder in structures of particle 
or molecules, e.g. the polydispersity of colloid particle 
sizes and the disorder in the stereoregularity, tacticity, 
and chemical structures of polymers. For case (i) the 
strength of disorder and random fields (S) is strongly de- 
pendent upon T, while for case (ii) the strength of the 
disorder is independent of T. The temperature-cooling 
paths drawn in figure [j] correspond to case (i) . 

Here we consider the non-Arrhenius behavior below T* n 
in more detail. We assume the divergence of the barrier 
height Eb arr i er between 'islands' is given by Eb arr i er = 
BT/(T-T ). This gives us the well-known Vogel-Fulcher 
(VF) law. Here we do not repeat the argument de- 
riving this relation (see e.g. §,0,0||J||l on the de- 
tails). Instead, we stress that the effective barrier domi- 
nating the structural relaxation E^ rier should be given 
by (i) EfJ Her = AG above T* and (ii) E e £ rier = 
AG + BT/(T - T ) below T£. This can be expressed 
analytically by the following modified VF law: 

r a =TZ>exp(AG/k B T)exp\f(T)B/k B (T-T )]. (4) 

Here f(T) is the Fermi-distribution-like function: f(T) = 
l/(exp(7(T — T^)) + 1), where 7 is a positive constant 
that controls the sharpness of the transition at T^. This 
modified VF law can naturally describe the crossover 
from the Arrhenius to non-Arrhenius behavior at 
and the divergence at To. The fitting of the tradi- 
tional VF law to data in a low temperature range of- 
ten produces an unrealistic attempt frequency (see 
e.g. p8[). This problem can be removed by using equa- 
tion (j4|) . From the condition of the disappearance of the 
non-Arrhenius behavior above T^, further, we obtain the 
relation B ~ AG(T* n - T a )/T* n . The fragility parameter 
[D] estimated by this relation with [D] = B/Tq is well 
correlated with D determined experimentally, as shown 
in the inset of figure 0. This provides us with a clear 
physical reasoning on the Angell plot JH|, including the 
correlation between bond strength S and the strong na- 
ture of liquids. 

Finally, we check the other main predictions of our 
model: (i) there occurs a dynamic transition from the 
non-cooperative (Arrhenius-type) to cooperative regime 
(Vogel-Fulcher- type behavior) at (near T m ); (ii) the 
disorder strength S, which can be correlated with the 
activation energy AG of the a relaxation above T^, de- 
termines the temperature distance between Tj^ and To; 

(iii) a stronger glass suffers from stronger disorder effects; 

(iv) an increase in pressure increases p, but decreases S 
for Ci < 0, simply because a locally favored structure 



having a greater volume is destroyed by applying pres- 
sure. This weakens the random disorder effects, which 
makes the system more fragile. Prediction (i) is quite 
natural and has already been recognized by many re- 
searchers (see e.g. [^8|), although there has so far been 
no physical reasoning on it. We stress that none of previ- 
ous models can explain this fact, since no models of glass 
transition have so far focused on crystallization phenom- 
ena themselves. This prediction is quite specific in the 
sense that T^ is a measurable quantity, while T c in MCT 
and T* in the theory of Kivelson et al. are not. Predic- 
tion (ii), on the other hand, can be confirmed by figure 
||, which indicates the strong correlation between AG (or 
AV) and T m — Tq. Prediction (iii) is consistent with : 
(a) the fact that the change in specific heat across T g is 
weaker for a stronger glass (see e.g. figure 4 of |Q); (b) 
the recent finding of Sokolov et al. [^9| on the stronger 
phonon scattering in stronger liquids; (c) a larger dis- 
tance between T g and To; and (d) the resulting larger (3k 
at T g for a stronger glass (see e.g. ||]). Prediction (iv) 
is also quite consistent with the experimental findings by 
Cook et al. |l],|30| that strong liquids become more frag- 
ile under a high pressure. This cannot be explained by 
the theory like MCT that is described by a single order 
parameter (density), since a plot of viscosity versus den- 
sity should be universal in such models. The relevance of 
our two-order-parameter (spin-glass) model to the glass 
transition is strongly supported by the validity of these 
predictions and also by the experimental results suggest- 
ing the similarity between structural glass and spin glass 

It is worth noting here the relation between our model 
and the model of Kivelson et al. |[l],[l2) , since both put a 
focus on 'frustration'. The most crucial difference is that 
we take a crystalline state (the ordering temperature T^) 
as the reference state, while they take a 'postulated' qua- 
sicrystal of a locally favored structure (the ordering tem- 
perature T*), which is prohibited in a real system due to 
frustration. In other words, they presuppose the avoid- 
ance of crystallization by hand. This leads to essential 
differences in physics, although both temperatures are 
claimed to be located near T m . The correlation of T^ to 
T m is quite natural for our case, while there seems to be 
no such justification for their case since we cannot expect 
any correlation between the hidden ordering point of a 
frustrated quasicrystal and that of a real crystal, which 
have essentially different symmetries. 

In summary, we propose a simple physical picture of 
glass transition on the basis of 'the two-order-parameter 
description of liquids' [B2], which connects structural 
glass and spin glass in a natural way. In our view, vit- 
rification is a result of the competition between density 
ordering and hidden local bond ordering. Our study in- 
dicates that frustrated systems such as spin glass and 
structural glass have the universal phase and dynamic 
behavior characterized by a complex energy landscape 
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peculiar to 'Griffiths-like' and 'spin-glass' phases. Here 
it should be noted that the fragile limit (S — > 0) with 
C2 = of the dynamic version of our model is mathemat- 
ically identical to the pioneering theory of Kirkpatrick 
and Thirumalai p0[ , which further has dynamic features 
similar to MCT [M (see the discussions in @). Thus, 
our model may give a clue to solve the problem of the 
existence of hidden disorder in mode-coupling equations 
pointed out by Bouchaud et al. Finally we stress 
that our model can provide us with a universal descrip- 
tion of the glass transition covering its strong to fragile 
limit. The strong nature of liquids increases with an in- 
crease in the disorder strength against density ordering 
(ciS and C2S). 
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FIG. 1. Schematic phase diagram of the liquid-glass 
transition. S is a measure of disorder strength against 
density ordering. SG stands for spin-glass phase. Crys- 
tallization occurs around T^, if it is not kinetically 
avoided. 
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FIG. 2. Relation between AG and T m - T . The data 
are taken from [12]. 1: n-butyl benzene; 2: isopropyl ben- 
zene; 3: propyl carbonate; 4: n-propanol; 5: o-terphenyl; 
6: salol; 7: dibutyl phtalate, 8: s-triaphthyl benzene; 9: 
glycerol; 10: a-phenyl-cresol; 11: boron oxide. The curve 
is to guide the eye. The inset shows the relation between 
[D] and D (see the text on their definitions) . 
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